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Abstract 

In this work, we calculate the CP-averaged branching ratios and the polarization fractions of 
the charmless hadronic Be — ?• A2A3 decays within the framework of perturbative QCD(pQCD) 
approach, where A is either a hght ^Pi or ^Pi axial- vector meson. These thirty two decay 
modes can occur through the annihilation topology only. Based on the perturbative calculations 
and phenomenological analysis, we find the following results: (a) the branching ratios of the 
considered thirty two — )• A2A3, decays are in the range of 10~^ to 10~^; (b) Be — )• oi^i, 
KiK+ and some other decays have sizable branching ratios and can be measured at the LHC 
experiments; (c) the branching ratios of Be — )• A2(^Pi)^3(^Pi) decays are generally much larger 
than those of Be — )• ^2(^Pi)^3('^Pi) decays with a factor around (10 ~ 100); (d) the branching 
ratios of Be — ?• KiK^ decays are sensitive to the value of 9k, which will be tested by the running 
LHC and forthcoming SuperB experiments; (e) the large longitudinal polarization contributions 
govern most considered decays and play the dominant role. 
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I. INTRODUCTION 



From the point of structure, the meson is a ground state of be system: which is 
hkely an intermediate state of the cc and 66-quarkonia, but should be very different from 
both of them since Be meson carries flavor B = —C = 1. When compared with the 
heavy-light Bg meson with q = {u,d,s), on the other hand, the decays of the B^ meson 
must be rather different from those Bu/B^/ Bg mesons since here both h and c can decay 
while the other serves as a spectator, or annihilating into pairs of leptons or light mesons 
(such as i^'+vr", etc). Physicists therefore believe that the B^ physics must be very rich 
if the statistics reaches high level [ll-[3|. In recent years, many theoretical studies on 
the production and decays of B^ meson have been done [2I, [sj, based on for example 
the Operator Production Expansion ||, NRQCDtsj], QCD Sum RulesQ, SU(3) flavor 
symmetry [3], ISGW II modeljsj, QCD factorization approach [9|, and the perturbative 
QCD (pQCD) factorization approach [lol- 13 1. 

On the experimental side, it is impossible to find a pair of B^B^ in the B-factory 
experiments (BaBar and Belle) since its mass is well above 6 GeV. Although the first 
observation of approximately 20 Be events in the Be — ?■ J /'^Iv decay mode was reported in 
1998 by the CDF collaboration f?], it was not until 2008 that two confirming observations 
in excess of So" significance were made by CDF and DO collaboration jl3] at Tevatron via 
two decay channels: the hadronic B^ — )■ J /'^ir^ decay and the semileptonic Be — )■ J /'^l^vi 
decay. 

At the LHC experiment, specifically the LHCb, one could expect around 5 x 10^" Be 
events per year0,l3|. And therefore, besides the charmed decays with large branching 
ratios, many rare Be decays with a decay rate at the level of 10~^ to 10~^ can also 
be measured with a good precision at the LHC experiments [7]. This means that, many 
Be — )■ hih2 decays ( hi are the light scalar(S), pseudo-scalar (P), vector(V), axial- vector (A) 
and tensor(T) mesons, made of light m, rf, s quarks ) can be observed experimentally. In 
the SM, such decays can only occur via the annihilation type diagrams. The studies on 
these pure annihilation Be decays may open a new window to understand the annihilation 
mechanism in B physics, an important but very difficult problem to be resolved. 

In 2004, by employing the low energy effective Hamiltonian |15,] and the pQCD ap- 
proach [iMi, we studied the pure annihilation decays B^ — )■ vrvr and presented the pQCD 
prediction for its branching ratiof]^: BrjBs — )• 7r"'"7r~) = (4.2 ± 0.6) x 10~^, which was 
confirmed by a later theoretical calculation [20| and by a very recent CDF measurement 



with a significance of 3.7a[2ll: Br{Bs tt+tt-) = (5.7 ± 1.5 ± 1.0) x This good 

agreement encourage us to extend our work to the case of Be decays. Although the charm 
quark c is massive (relative to the known light quarks u, rf, and s), the Be meson has 
been treated as a heavy-light structure in this work because of the ratio me/mB^ ~ 0.2, 
which means that the large part of the energy is carried by the much heavier b quark in 
a Be meson. With this assumption, we also employ the fcy factorization theorem to the b 
decay in Be meson, in a similar way as for the decays of B^ and B^ mesons. 

During past two years, based on the pQCD factorization approach, we have made a 



systematic study on the two-body charmless hadronic decays of Be — )■ PP, PV,VV [10 
Be ^ SP, SV [nj and Be AP, AV For all the considered pure annihilation Be 

decay channels, we calculated their CP-averaged branching ratios and longitudinal polar- 
ization fractions, and found some interesting results to be tested by the LHC experiments. 
In this paper, we extend our previous investigation further to the charmless hadronic 
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Be — AA decays. The axial-vector mesons involved are the following: 

ai(1260), 6i(1235), fs:i(1270), 7^1(1400), /i(1285), /i(1420), /ii(1170), /ii(1380). (1) 

All the thirty two decay modes are the pure annihilation decay processes in the SM. 

The internal structure of the axial-vector mesons has been one of the hot topics in 
recent years |22N24|. Although many efforts on both theoretical and experimental aspects 



have been made [25H30| , we currently still know little about the nature of the axial-vector 
mesons. Our study will be helpful to understand the structure of these mesons. 

As one of the popular factorization tools based on the QCD dynamics, the pQCD 
approach can be used to analytically calculate the annihilation type diagrams. Besides 
the good agreement between the pQCD prediction and the newest CDF measurement for 
Br{Bs TT+TT-), the pQCD prediction of Br{B^ D-K+) ^ (4.6 ± 1.0) x IQ-^ for 
the pure annihilation B^ decay as presented in Ref. 3l|] also be consistent well with the 
data |30|. We therefore beheved that the pQCD factorization approach is a powerful and 
consistent framework to perform the calculation for the annihilation type By^^d,s decays, 
and extend our work to the cases of B^. decays. 

The paper is organized as follows. In Sec. [TTl we give a brief review about the axial- 
vector meson spectroscopy, and the theoretical framework of the pQCD factorization 
approach. We perform the perturbative calculations for considered decay channels in 
Sec. mil The analytic expressions of the decay amplitudes for all thirty two Be — )■ AA 
decays are also collected in this section. The numerical results and phenomenological 
analysis are given in Sec. IIVI The main conclusions and a short summary are presented 
in the last section. 



II. THEORETICAL FRAMEWORK 
A. Axial-vector mesons and mixings 

In the quark model, there exist two distinct types of light parity-even p-wave axial- 
vector mesons, namely, ^Pi{Jpc = l^"*") and ^Pi{Jpc = l"*"") states: 

nonet :ai(1260),/i(1285),/i(1420) and Kia; 

^Pi nonet :6i (1235), /ii (1170), /ii (1380) and K^b- (2) 

In the SU(3) flavor limit, the above mesons can not mix with each other. Because the 
s quark is heavier than m, d quarks, the physical mass eigenstates -ft'i(1270) and -ft'i(1400) 
are not purely ^Pi or ^Pi states, but believed to be mixtures of Kia and -ft'ijll- Analogous 
to T] and rj' system, the fiavor-singlet and flavor-octet axial-vector meson can also mix 
with each other. 

The physical states Ari(1270) and ^1(1400) can be written as the mixtures of the Kia 
and KiB states: 

^^1(1270) = sitiOkKia + cos OkKib , 

Ki{im) = cosOkKia - sinOKKiB , (3) 

^ For the sake of simplicity, we will adopt the forms ai, bi, K , K , f , f , h and h to denote the 
axial-vector mesons ai (1260), 6i(1235), /s:i(1270), Xi(1400), /i(1285), /i(1420), /(,i(1170) and /ii(1380) 
correspondingly in the following sections, unless otherwise stated. We will also use Ki, f[ and h[ to 
denote fCi(1270) and is:i(1400), /i(1285) and /i(1420), and /ii(1170) and /ii(1380) for convenience 
unless otherwise stated explicitly. 
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where 6k is the mixing angle to be determined by the experiments. But we currently have 
little knowledge about 9^ due to the absence of the relevant data, although it has been 
studied for a long time |22|-|24| In this paper, for simplicity, we will adopt two reference 
values as those used in Ref. [2J]: 9k = ±45°. 

Analogous to the rj-rj' mixing in the pseudoscalar sector, the /ii(1170) and /ii(1380) 
(^Pi states) system can be mixed in terms of the pure singlet hi and octet hg, 



/ii(1170) = sin 6*1 hs + cos 6*1 hi, 
/ii(1380) = cos^^i hs — sin6'i hi. 

Likewise, /i(1285) and /i(1420) (the ^Pi states) will mix in the same way: 

/i (1285) = sin ^3/8 + cos ^3/1, 
/i(1420) = cos^3/8-sin^3/i. 

where the flavor contents of hi^g aiid /i,8 can be written as 

hi = fi = —= [uu + dd + ss) , 

hs = fs = —j= (uu + dd — 2ss) . 
v6 

The values of the mixing angles ^1^3 can be chosen as j24| : 

9i = 10° or 45°; 9^ = 38° or 50°. 



(4) 



(5) 



(6) 



(7) 



B. Formalism 



In the pQCD factorization approach, the four annihilation Feynman diagrams for Be — )■ 
A2A3 decays are shown in Fig.l, where (a) and (b) are factorizable diagrams, while (c) 
and (d) are the non-factorizable ones. The initial b and c quarks annihilate into u and d/s, 
and then form a pair of light mesons by hadronizing with another pair of qq {q = {u, d, s)) 
produced perturbatively through the one-gluon exchange mechanism. Besides the short- 
distance contributions based on one-gluon-exchange, the qq pair can also be produced 
through strong interaction in non-perturbative regime (final state interaction(FSI), for 
example) . 



B 




(a) (b) (c) id) 



FIG. 1. The annihilation Feynman diagrams for Be — >■ A2A3 decays. (a) and (b) are factorizable 
diagrams; while (c) and (d) are the non-factorizable ones. 
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For the considered — )■ A2A3 decays, the key point is to calculate the corresponding 
matrix elements: 

M oc< A2A3\neff\B^ > (8) 
where the weak effective Hamiltonian "He// is given by [l^ 

with the current- current operators O12, 

O2 = U^ril - 75)i^/3C,7'^(l - 75)&a , (10) 

where Vcb, Vud (D = d, s) are the CKM matrix elements, Cj(/i) are Wilson coefficients at 
the renormalization scale fi. 

Although the dominance of the one-gluon exchange diagram seems favored by the data 
of — TT~^Tr~ and 5° — )■ D~K^ decays, according to the good a gree ment between our 
calculations based on the pQCD approach [l9, 31| and the data [21, sO], we currently still 



do not know whether the short-distance or the non-perturbative contribution dominate 
for Bf. annihilation decays. We here first assume that the short-distance contribution 
is dominant, and then calculate the matrix element in Eq. ([8]) by employing the pQCD 
approach, provide the pQCD predictions for the branching ratios and longitudinal polar- 
ization fractions, and finally wait for the test by the LHC experiments. 

We work in the frame with the B^. meson at rest, i.e., with the B^. meson momentum 
Pi = -^(1, 1, Ot) in the light-cone coordinates. We assume that the A2 (A3) meson moves 
in the plus (minus) z direction carrying the momentum P2 (P3) and the polarization vector 
€2 (€3). Then the two final state meson momenta can be written as 

P. = ^{l-rlrlOr), P, = '^{rll - rl^). (H) 

where r2 = m^j/m^, and = mA-i/friB- The longitudinal polarization vectors, €3 and 
63 , can be defined as 

4 = ^^(1-4-40^), e',=^^{-rll-rlOr). (12) 
v2mA2 v2mA3 

The transverse ones are parameterized as = (0, 0, It), and = (0, 0, 1^). Putting the 
(light-) quark momenta in B^, A2 and ^3 mesons as fci, k2, and k^, respectively, we can 
choose 

fci = (a;iP+,0,kiT), k2 = {x2P^,0,k2T), k^ = {0,XsP,- ,ksT). (13) 
Then the decay amplitude can be written conceptually as the following form, 



M{B,^ A2A3) = < A2A3\neff\B,> jc 



dx 1 dx2 dx^hi dhi 62 dh2 db; 



'3 



X 



Tr [C(t)$B,(xi, 61)$^, (X2, &2)$A3 {X3, bs)H{x,, 6„ t)St{x,) e'^^*)] (14) 
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where bi is the conjugate space coordinate of kix, and t is the largest energy scale in 
function H{xi, bi,t). The large logarithms ln{mw /t) are included in the Wilson coefficients 



C{t). The large double logarithms (In Xj) are summed by the threshold resummation [32 



and they lead to St{xi) which smears the end-point singularities on Xi. The last term, 
g-5(i)^ is the Sudakov form factor which suppresses the soft dynamics effectively [s^. Thus 
it makes the perturbative calculation of the hard part H applicable at intermediate scale, 
i.e., thb^ scale. We will calculate analytically the function H{xi,bi,t) for the considered 
decays at leading order(LO) in expansion and give the convoluted amplitudes in next 
section. 



III. THE DECAY AMPLITUDES IN THE PQCD APPROACH 

For an axial-vector meson, there are three kinds of polarizations, namely, longitudinal 
(L), normal (A^), and transverse (T). The Be — )■ ^2(62, -P2)^3(e3, -P3) decays are charac- 
terized by the polarization states of these axial-vector mesons. 



A. Decay Amplitudes with different polarization 

The decay amplitudes J^h are classified accordingly, with H = L,N, T, 

Mh = (ml Ml, mlMNel{T)eliT), iMTe''^'^''el^{T)elp{T)P2,P,p) . (15) 

where e(T) stands for the transverse polarization vector and we have adopted the nota- 
tion e*^^^^ = 1. Based on the Feynman diagrams shown in Fig. [H we can combine all 
contributions to these considered decays and obtain the general expression of total decay 
amplitude as follows, 

MniB, ^ A2/I3) = V:,VuD {fB.Ff^'^'a, + M^l^^^^Ci} , (16) 

where ai = Ci/3 + C2 1, while F^aSi' ^na^ dcuotc the Feynman amplitudes with 
three polarizations for factorizable and nonfactorizable annihilation contributions, respec- 
tively. 

The explicit expressions of the function F^^.^ and M^^f^ in the pQCD approach can 
be written as the following form: 

nl POO 

Fj^ = SnCprn^^^ / dx2dxs / b2db2b^db^ 

X { [a;202(x2)03(a;3) + 2r2r3 ((sa + 1)02(3^2) + {^2 - l)(t)\{x2)) (t)i{x^)\ 

xEfa(ta)hfa{l - X3,X2, 63,^2) + E fa(tb)h fa{x2, 1 - X^MM) 

X [(X3 - I)(p2{x2)(pz{x^) + 2r2r30^(x2) ((X3 - 2)0^(X3) - X30^(X3))] } , (17) 



^ One should note that ai here just stands for the combined Wilson coefBcient, not the abbreviation for 
axial- vector meson ai(1260). 
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-J "Jo Jo 

X {[{Tc -X3 + 1 ) 02 (a;2) 03(3:3) + '^2'^3 (02(2:2) ((3rc + X2- X3 + 1) 
X0^(X3) - (rc - X2 - X3 + 1)03(X3)) + 02(x2)((rc - X2 - X3 + 1)03(X3) 
+ {rc -X2 + X3- 1)03(X3)))] Ena{tc)Ka{^2, 61, 62) 

- [{rb + rc + X2- 1)02 (3:2) 03(3:3) + r2r3 (0^(0:2) ((4rb + + X2 - 0:3 - 1) 

X03(X3) - (rc + X2 + X3 - 1)03 (X3)) + 02(x2)((rc + X2 + X3 - 1)03(^3) 

-{rc + X2-X3-l)(p\{x3)))]Ena{td)h'^nai.X2.XzMM)] , (18) 

/I /'OO 
dx2dx3 / b2db2b3db3 r2r3 

X {[(0:2 + 1)(0^(X2)0^(X3) + 0^(X2)0^(X3)) + (0:2 - l)(0^(a;2)0^(a;3) + 02(^2)0:^(x3))] 

xEfaita)hfa{l - X3, X2, 63, 
+ [(X3 - 2)(0^(a;2)0^(x3) + 0^(X2)0;;(X3)) - X3 (0^(X2)0^(X3) + 0^(X2)0^(X3))] 

xEfa{tb)hfa{x2, 1 - X3, 62, ^3)} , (19) 

^na= o ^C'irml^ / dX2dX3 / 61^6162(^62 '^2?'3 

<J "Jo Jo 

X {re [0^(X2)03(X3) +0^(x2)03(x3)]E„„(tJ/l^„(x2,X3,6i,62) 

-rfc[0»(a;2)0^(x3) + 0^(x2)0^(x3)]i^;„a(td)/^i(a:2,X3,6i,62)} , (20) 



FJ^ = IQ-KCpm^^ I dx2dx3 I 62^6263^63 r2r3 
'0 Jo 



f 1 ^00 

7. ' 

X {[(X2 + 1)(0^(X2)0:^(X3) + 0^(X2)0^(X3)) + (^2 - 1) (0^(x2)0^(x3) + 0^(X2)0:^(X3))] 

xEfa{ta)hfa{l - X3, X2, 63, 62) 
+ [{X3 - 2)(0^(a;2)0^(x3) + 0^(X2)0^(X3)) - X3 mX2m^3) + <p2i^2)4>lix3))] 

xEfa{tb)hfa{x2, 1 - X3, 62, 63)} , (21) 

Mna = ^^T^Cprnj^l dx2dx3 / 616/6162^62 r2r3 

X {re [02(2:2)03(2:3) + 02(2:2)03(2:3)] Ka(tc)/iL(a:2, 2:3, &1,&2) 

-rfe[0^(x2)0Ka;3) + 0^(x2)0^(a;3)]i?„,a(td)/ii(2:2,X3,6i,62)} . (22) 

where r2(3) = mA^^^^/^Ba^ '^c(b) = ^c{b)/TTT'Bc- The exphcit expressions for the distribution 
amphtudes 0^, 0^, 0^, 0^, 0^ and 0^ are given in the Appendix A. The definitions and 
expressions of the hard functions (hfn, hnn), {Efa,Ena) and hard scales (ta,tb,tc,td) can 



be found in Appendix B of Ref. [10| and references therein 



B. Decay Amplitudes for the considered decay modes 

Now we can write down the total decay amplitudes for all thirty two Be — )■ A2A3 
decays. The decay amplitudes of the sixteen AS* = decay modes are the following: 

V2Mh{B. ^ a+a?) = KtK. {fs. [pfj- - F;>^) a, + [mS^ - M^^^) Ci}(,23) 
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v/2M,(S. ^ btbl) = K;K. {fB. - F«) a, + (M:i;> - M«) C,} ,(24) 

V2Mh{B. ^ a+6?) = KlKd {/i.. (F^i;^!" - F;!;'^^) a^ + (mI^J" - M^^^^^) C,} ,(25) 
v^M^(Se ^ 6+a?) = KtKd {/b. - F<^) a, + (M:i< - M<;'^) C,} (26) 



sin 6'3 



(27) 



sin 

71~ 



/Be ( Ffj + F|;i) ai + [Mis + M£;i) 



COS ^3 



(28) 



Mh{B, ^ btf') 
Mh{B, -> 6+/") 
Mh{Bc a^h') 
Mh{Bc a+h") 
MniB^-^hth') 
MHiB.^bth") 



- Mh{Bc atf'){ai bi), 
: Mh{Bc atf"){ai 6i), 
■.MH{Bc^atf'){f ^h,es^ei), 

■.MH{B,^atf"){f ^h,es^e,), 

■■ Mh{Bc ath'){ai bi), 
■ Mh{Bc ath"){ai 6i), 



(29) 
(30) 
(31) 



Mh{B. ^ K'"K' ■ ) = V^V^, [sin' 9^ {fsP'^f'^a, + M^f^^C,) 

- cos Ok sin Ok [fsP'^f'-a, + M^-^-C,) 
+ cos Ok sin (/B.Fj-^-ai + M^^f^^C,) 
cos' Ok [Ib^f'^^h 



«1 + ^^^^na;i? 



(32) 



M^,(S, ^ i^' V"+) = KIK. {-COS Sin Ok {fB.F^f'^a, + M^^f^^C,) 



+ sin^ f /B.F.^^-^^^ai + Mf -^-Ci 



+ cos^ Ok {fBFf^f'^a, + M^f^^C,) 
- cos sin Ok (fs^F'^lf'-a^ + M^f^^C, 



(33) 
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Mh{B. ^ X^V'+) = KtKd {cos sin {fsF^f'^a, + M^f^^C,) 

+ cos Ok sin Oj, [fsFj^^f'^a, + M^^-^i) } , 



- COS sin Ok {fsF^f'^a, + M^f^^C,) 
+ cos^;,sin^^ [fBF'^,fa, + M^f^-Ci) 

- sin^ 9k {fB^F^f'^a, + M^f^'^C^) } ; 

The decay amplitudes of the sixteen A^* = 1 decay modes are of the form: 



Mh{B, ^ K' at) = V2Mh{B, K a'i) 



^cosOk [fBFf}-/a, + M^!.fc,]] , 
- sin eK[fB.Ff!fa, + Mf^.fc,] } , 



Mh{B^ K'%i) = V2Mh{Bc K'^bl) = Mh{Bc ^ ir'V)(«i ^ &i) , 



X 



sin ^3 sin 9 k 
cos 9s cos 9 k 
cos 9k sin 6*3 



f ( TT^'^B, 
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Mh{Bc ^ K'^f") = KtK. 



sin 03 sin r f K^Afr , p/f^iA\ ^ , f m^^^^^ ^ M^^^'A r 

[IB, [^fa-H + ^fa;H j «1 + [^^^na;H + ^^^na;H ) ^1 

COS 93 sin Ok U /pKiA/g" op/I^iaA , /-^-^ia/s" r,^fiKrA\^' 



r cos 93 cos 9 k r / Ki^/r /f \n 
^ I ^ [jBc [^fa;H + ^ fa;H j «1 + (^^^na;// 

I cos^j^sin^s r f x.^fS opfi^^A\ ^ , f j^K.^fS c.^. 



na;H + ^^-^na;H J '-'1 



V6 

COS 93 sin 9k {, /j^. 
r ^ 



J Be [^fa;H 



sin 9k sin ^3 r / i^isj 
V6 L-^""^ V ^"'"^ 



yjB^y- fa-^H +^fa;H ) + [^na;H + ^^na;H j 

fa;H - ^^fa;H J «1 + y^^na;H " ^^^^na;H ) 



I COS ^3 COS ^K r. /pXiA/g" Qpfi^^A^ , //r^-^iA/g" 9 ;j^/|i^lAX^ i 
UBA^fa;H - fa;H ^ + V^^na;H - ^^^^na;H J<^lJ 

I sin^ssin^j^ /pi^is/r , pfl^is^ , / a^^is/i" , M^^^^''\r 1 

^^^1 UBA^fa;H + ^ fa;H J«l + l^^^na;if + ^^^na;H J<^lJ 

COS ^3 sin .pi^is/g" npfil<^B^ ,(.rK,Bft 
UBA^fa;H ' fa;H ^ + V^La;H 



Mh{Bc ^ K"'h') = Mh{Bc ^ ^ /i, ^3 ^ ^i) , 

MniBc ^ iiT'^/i") = Mh{Bc ^ ^ /i, ^3 ^ ^1) , 

MniBc ^ = A4if(5c ^ ^ /i, ^3 ^ ^1) , 

Mh{Bc ^ = A^i/(5e ^ K"^f"){f ^h,93^ 9i) . 



, (42) 



■2mS-)c4, (43) 



(44) 
(45) 
(46) 

(47) 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we will calculate numerically the BRs and polarization fractions for 
those considered thirty two Be — )> A2A3 decay modes. First of all, the central values of 
the input parameters to be used are the following. 
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Masses (GeV): 



mw = 80.41; m^^ 
= 1.23; uik^^ 
irtb^ = 1.21; mKj_B 



Decay constants (GeV): 

fa, = 0.238 

fb, = 0.180 
/b. = 0.489 



6.286; 

1.32; 

1.34; 



0.250; 



fK,, = 0.190; 



rrib 



fhi 



4.8; 

1.28; 

1.23; 



0.245; 
0.180; 



rric = 



fh ■ 

fhs 



1.5; 
: 1.28; 
= 1.37; 



(48) 



0.239; 
: 0.190; 



QCD scale and Be meson lifetime: 



A. 



(^4) 
'MS 



0.250 GeV, tb, = 0.46 ps. 



(49) 



(50) 



For the CKM matrix elements we use A = 0.814 and A = 0.2257, p = 0.135 and f] = 0.349 



34( 1 . In numerical calculations, central values of input parameters will be used implicitly 
unless otherwise stated. 

For these considered Be — )■ A2A3 decays, the decay rate can be written explicitly as. 



r 



167rm^ 



(51) 



=L,T 



where |Pc| = |P2z| = iPszl is the momentum of either of the outgoing axial- vector mesons 



The polarization fractions /l(||,±) can be defined as [35 

where the amplitudes Ai{i = L, ||, ±) are defined as, 

Al = -^ml^Mi, All = ^V2mlMN, A± = ^mA^mA^VW^^Mr , (53) 
for the longitudinal, parallel, and perpendicular polarizations, respectively, with the nor- 

Ps/imA^ rriAs). These 



malization factor ^ = y G|,Pc/(16'7rm^^r) and the ratio r = P2 
amplitudes satisfy the relation. 



\Al\' + \Aii\' + \A^\' 



1 . 



(54) 



following the summation in Eq. (15T|) . 

By using the analytic expressions for the complete decay amplitudes and the input 
parameters as given explicitly in Eqs. (123|) - (!50|) . we calculate and then present the pQCD 
predictions for the CP-averaged BRs and longitudinal polarization fractions (LPFs) of 
the considered decays with errors in Tables HHVl The dominant errors arise from the 
uncertainties of charm quark mass rric = 1.5 ± 0.15 GeV and the combined Gegenbauer 
moments of the axial-vector meson distribution amplitudes, respectively. 
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TABLE 1. The pQCD predictions of BRs and LPFs for Be 
of the dominant errors is explained in the text. 



(ai, 6i)(ai, 61) decays. The source 



AS = 
Decay modes 


BRs (10-^) LPFs (%) 


AS = 
Decay modes 


BRs (10-5) LPFs (%) 


Be — ?" a^a^ 
Be ^ a+b? 


0.0 


Be ^ b+b? 

Be hf4 


0.0 



A. The pQCD predictions for AS = decays 

In Table [T] and II, we show the pQCD predictions for the branching ratios and the 
longitudinal polarization fractions of the sixteen AS = decays. 

For both the — t- afa^ and 6^65 decays, since the quark structure of and 
are the same one , {uu — dd) / a/2, the contributions from uu and dd components to the 
corresponding decay amplitude as shown in Eqs.f l23|l2^ will interfere destructively, and 



therefore will cancel each other exactly at leading order and result in the zero BRs for 
these two channels, as illustrated in the Table I. For the possible high order contributions, 
they will also cancel each other due to the isospin symmetry between u and d quarks. As 
for the non-perturbative part, we currently do not know how to calculate it reliably. But 
we generally believe that it is small in magnitude for B meson decays. Consequently, we 
think that a nonzero measurement for the branching ratios of these two decays may be a 
signal of the effects of new physics beyond the SM. 

For Be — )■ a^bi and Be — )■ b^a^ decays, however, the pQCD predictions for their BRs 
are rather large, as given in Table [T] 



Br{Be -> atb\) = Br{Bc b^a^) ^ 2.2 x 10" 



(55) 



Besides Be — ?► a^b^ and 6]^a° decays, other six = decays, such as the Be — )■ 
and Be — )■ K^K^ decays, also have a large branching ratios at the 10"^ level, as listed 
in Table II. According to the studies in Ref. these Be decay modes with a branching 
ratio at 10"^ level could be measured at the LHC experiments 0]. 

Besides the large branching ratio at 10"^ level, the Be — t- K^Kf decay modes also 
have a strong dependence on the value of the mixing angle Ok-, as shown by the numbers 
in Table II. If these channels are measured at LHC experiments with enough precision, 
one can determine the 9k by compare the pQCD predictions with the data. In order to 
reduce the effects of the choice of input parameters, we define the ratio of the branching 
ratios between relevant decay modes: 



Br{Be ^ iri(1270)°iri(1400)+)pQCD 
Br{Be ^ :^i(1270)o/^i(1270)+)pQCD 



3.0, for Ok = 45° , 
0.7, for Ok = -45° ; 



(56) 



BrjBe -> iri(1270)°iri(1400)+)pQCD 
Br {Be ^ (1400)07^1 (1270)+)pQCD 



2.0, for Ok 
0.5, for 9k 



45° , 
-45° ; 



(57) 
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TABLE II. Same as TableHbut for (a+, h[) decays. 



AS = 
Decay modes 


fl., — 38° 

OO 

BRs (10-6) LPFs (%) 


— 50° 

BRs (10-6) LPFs (%) 


ai(1260)+fi(1285) 
Be ai(1260)+fi(1420) x IQf. 


6.5+?,S(m^)+?n(a,) 83.6+^-f 


6.1+iSfm^)+no(a,) 84.0+H 

— 0.9V c/_o.9v */ —4.0 

3.9t°:^(m,)ll:3(a,) 78.5ll3^9 


AS = 
Decay modes 


BRs (10-^) LPFs (%) 


OKJ 

BRs (10-^) LPFs (%) 


bi(1235)+fi(1285) 
Be ^ bi(1235)+fi(1420) 


2.8+^-km,)+i-«(a,) 65.2+?^-3 

—0.5 V — 0,9 v^/ —16.4 
1.4t°:2(me)l°:^(a,) 100.0 ±0.0 


3.0+^-^(m^)+i-5(a,) 68.7+?j,-I 

— 0.8V cy_o,9v^/ ' —14.6 

1.2tli{m,)lli{a,) 100.0^°:° 


AS = 
Decay modes 


9i = 10° 
BRs (10-6) LPFs (%) 


6*1 = 45° 
BRs (10-6) LPFs (%) 


B^ ai(1260)+hi(1170) 

Be ai(1260)+hi(1380) x 10 


1.3+^-^('m^)+^-5('a,) 86.3+q2 

O.oV c; — 0.4V ""^—9.8 


0.7+^-^(m^)+^-?(a,) 73.1+L'ic 

" ' — 0.4v "^c/— 0.2^""*/ '"-'■-^ — 28.8 

6.8t°;2(mc)l|I(ai) 100.01?;° 


A <? — n 
Decay modes 


—10° 

(7]^ — lU 

BRs (10-5) LPFs (%) 


BRs (10-5) LPFs (%) 


Be bi(1235)+hi(1170) 
Be ^ bi(1235)+hi(1380) 


8.1+^i(me)+H(ai) 96.4+j-° 
2.5t[]:?(me)tJ:^(a,) lOO.OtO;" 


10.3+H(mc)+tI(ai) 96.4+?i 
0.3l[!:^(mc)l[]:^(ai) 100.0 ±0.0 


A5 = 
Decay modes 


= 45° 

BRs (10-5) LPFs (%) 


Ok = -45° 
BRs (10-5) LPFs (%) 


Be ^ Ki(1270)°Ki(1270)+ 
Be ^ Ki(1400)°Ki(1400)+ 


1.2l°;?(me)lJJ(a.) 99.7^?:^ 
2.8li:g(me)iy(a,) 72.7115.8 


2.9tl:^(m,)l|^(a,) 71.9ti6-2 
l.ltroMlliia,) 99.7t?;° 


Be ^ Ki(1270)0Ki(1400)+ 
Be ^ Ki(1400)°Ki(1270)+ 


3.7l}i(me)l|i(a,) 96.2l|5 
1.9t°:5(me)l?i(a,) 94.6l?o6^ 


1.9«J(m,)l?i(a,) 94.8+?o'3 
3.7li:?(m,)l|?(a,) 96.lt|6 



^ Here, the factor 10 is specifically used for the BRs. The following one has the same meaning. 

BrjB, ^ Zi(1400)°i^i(1400)+)pQCD _ f 2.3, for 9k = 45° , 
Br{B, ^ (1270)0/^1 (1270)+)pQCD I 0.4, for Ok = -45° ; 

The LHC experiments can measure these ratios with a better precision than that for a 
direct measurement of branching ratios for individual decays. We suggest such measure- 
ments as a way to determine the mixing angle 9k at LHC. 

B. The pQCD predictions for AS" = 1 decays 

In Table III, IV and V, we show the pQCD predictions for the branching ratios and 
the longitudinal polarization fractions of the sixteen AS* = 1 decays. 

First of all, when compared with those AS* = decays, these AS = 1 decays are CKM 
suppressed due to the factor ~ 0.04, as can be seen easily from the expressions 

for the decay amplitudes as given in Eqs. fl23l) to (HTll . The pQCD predictions for the 
branching ratios of these Be decays are at the level of 10-6 to 10-*^ , much smaller than 
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TABLE III. Same as Table H] but for Be — s- Kiai, Kibi decays. 



A5 = 1 
Decay modes 


Oj^ = 45° 
BRs (10-^) LPFs (%) 


= -45° 

BRs (10-^) LPFs (%) 


Ki(1270)°ai(1260) + 
Be Ki(1400)°ai(1260)+ 


4.6+}-^(m,)+H(a,) 79.2+}2-4 

— l.Ov c/— 2.4V"'«/ ' — 16.3 

8.0li-3(m,)+|^(ai) WO.Otll 


8.3+}-^(m.)+H(a,) 99.3+9-f 

— 1.8V c/— 3.9V"*/ '^'^ 5.5 

4.5ll:?(m,)iy(a,) 8l.3t\li 


Be ^ Ki(1270)+ai(1260)° 
Be Ki(1400)+ai(1260)° 


^■Ot',lMt\i{a,) 100.0^°;° 


^.2tft{m,)tli{a.) 99.3t0| 
2.2l°|(m,)l?:?(a.) 81.3l}2-5 


AS= 1 
Decay modes 


9k = 45° 
BRs {W'^) LPFs (%) 


9k = -45° 
BRs (lO"'') LPFs (%) 


Be ^ Ki(1270)0bi(1235)+ 
Be ^ Ki(1400)°bi(1235)+ 


l-6t'oiMlli{a.) 91.3+^;? 
1.3t°-^(me)l|]|(ai) 100.0 ±0.0 


1-<2K)T7(«0 100.0^°;° 
1.5l°|(m,)iJ|(a.) 93.6ii:? 


Be ^ Ki(1270)+bi(1235)° 
Be ^ Ki(1400)+bi(1235)0 


0.8tS:^(m.)l°|(a.) 91.411? 
0.7+0|(m,)l°:^(ai) 100.0 ± 0.0 


o.7^g:?K)l°i(a.) ioo.oi°:° 

0.8l0J(m,)10|(a,) 93.6i|] 



that for those AS = decays. Most of them, for example — )■ i^iOi and i^i/i decays 
with BRs around 10"^ or less, are hardly to be detected even at the LHC experiments. 

For the — ?• i^i&i decays, the pQCD predictions for the BRs are in the order of 10~^, 
much larger than the BRs of the B^ Kiai decays, since the ^Pi meson behaves very 
different from the ^Pi state. From the numerical values in Table III, we can also define 
the following ratio 

Br{Be Ki(1270)°6^") pQCD Br{Bc — )■ i^i(1400)°6+)pQCD r, 

Br{B, ^ iri(1270)+6?)pQCD 5r(5, ^ ifi(1400)+6?)pQCD ^ ' 

for both 9k = ±45°. Such decays have a weak dependence on the variation of Ok- 

In Table IV, we show the pQCD predictions for the BRs and LPFs for B^. — )■ K^f[ 
decays with 6*3 = 38° (1st entry) and 6*3 = 50° (2nd entry), respectively. In Table V, 
similarly, we show the pQCD predictions for the BRs and LPFs for Be — ?■ K^h[ decays 
with 61 = 10°(lst entry) and 63 = 45°(2nd entry), respectively. 

One can see from the numerical results in these two tables that all Be — ?■ Ki{f[, h[) 
decays have a weak or moderate dependence on the mixing angles 61 and 63. It is difficult 
to measure 61 and 63 through the considered Be decays. 

For Be — 7- K^hi{1380) decays, the pQCD predictions for their BRs show a relatively 
strong dependence on the mixing angle 6k- The LHC measurement of these decays may 
also help to constrain the size and sign of 6k- 

Frankly speaking, the theoretical predictions in the pQCD factorization approach still 
have large theoretical errors induced by the large uncertainties of many input parameters 
and the meson distribution amplitudes. Any progress in reducing the error of input 
parameters will help us to improve the precision of the pQCD predictions. 

It is worth of stressing that we here calculated only the short-distance contributions in 
the considered decay modes and do not consider the possible long-distance contributions, 
such as the rescattering effects, although they may be large and affect the theoretical 
predictions. Strictly speaking, it is the task after the first measurements of the B^ meson 
decays and thus beyond the scope of this work. 
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TABLE IV. Same as Table H but for K:^ f[ decays with 6*3 = 38° (1st entry) and 63 = 

50° (2nd entry), respectively. 



AS = 1 
Decay modes 


Ok = 45° 
BRs (10-^) LPFs (%) 


Ok = -45° 
BRs (10-7) LPFs (%) 


Be ^ Ki(1270)+fi(1285) 


l-4lHK)lg:?(a.) 65.ltlli 
1.7tli{m,)tliia,) 69.ll?2.i 


l.5t'oiimc)tliia.) 92.ll?38o 


Be ^ Ki(1400)+fi(1285) 


1.5tUimc)tliia.) 92.11^2^8 


l.7lliim^)troia.) 


Be ^ Ki(1270)+fi(1420) 


0.9«:^(m,)l°J(aO 81.6^^3.5 
0.6^[Ji(me)+°:^(a.) 78.51^^:? 


AA+li{m,)+l',{a,) 71.5tti 
4.4l°J(me)71:^(a.) 73.2+^;^ 


Be ^ Ki(1400)+fi(1420) 


4.3l°-^(m,)l}-?(a,) 71.9^^1 
4.4l0:|(me)t}:^(a,) 73.6^^1 


0.9lSj(me)i°|(a.) 81.9tlli 



TABLE V. Same as Table U but for K^h[ decays 61 = 10°(lst entry) and 6*1 = 45°(2nd 

entry), respectively. 



A5 = 1 
Decay modes 


Ok = 45° 
BRs (10-6) LPFs (%) 


Ok = -45° 
BRs (10-6) LPFs (%) 


Be ^ Ki(1270)+hi(1170) 


1.4l°|(me)iJ:|(a,) 94.5i|3 
0.6t'oiim.n'M) 87.9ttt, 


1.6i|]:I(me)ti:°(a.) 98.51°:^ 
0.2t'o'oi^c)tli{a.) 92.9tli, 


Be ^ Ki(1400)+hi(1170) 


1.6tU(mc)t\-'M) 98.5t°:^ 
0.2t'olMt',l{a,) 93.0tlf,s 


lAtli{m,)tl-l{a,) 9A.6+li 
0.5l°:^(m,)+0:6(a.) 88.1^11, 


Be Ki(1270)+hi(1380) 


0.9t'oi{mcm{a.) 98.5l?i 
1.8l°i(me)lJi(a.) 98.61°:? 


1.5l°:l(mc)l|]:?(ai) 89.6l|^ 
2.8li:i(me)ll:i(a,) 94.3^^:^ 


Be Ki(1400)+hi(1380) 


l■5tliim.n|ia^) 89.8tl-l 
2.8tli{m,)t\i{a^) 94.41^:? 


0.9tli{m,ni{a.) 98.51?:^ 
l-7ig|(-c)l^:^(a.) 98.6^°:^ 



V. SUMMARY 



In this paper, we studied the thirty two charmless hadronic A2A3 decays by em- 

ploying the pQCD factorization approach. These considered decay channels can only occur 
via the annihilation type diagrams in the SM. The pQCD predictions for the CP-averaged 
branching ratios and longitudinal polarization fractions are analyzed phenomenologically. 

From our perturbative evaluations and phenomenological analysis, we found the fol- 
lowing results: 

1. The branching ratios of the considered thirty two — )■ AA decays are in the range 
of 10"^ to 10"^; Be — 7- cti&i, K^Kf and some other decays have sizable branching 
ratios (~ 10"^) and can be measured at the LHC experiments; 
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2. The branching ratios of Be A2{^ Pi) A^l^ Pi) decays are generally much larger than 
those of Be — )■ A2{^Pi)A3{^Pi) decays with a factor around (10 ~ 100) because of 
the rather different QCD behavior between ^Pi and ^Pi states; 

3. For Be — >■ AA decays, the branching ratios of AS* = processes are generally much 
larger than those of AS* = 1 ones. Such differences are mainly induced by the CKM 
factors involved: Vud ~ 1 for the former decays while Vus ~ 0.22 for the latter ones. 

4. The branching ratios of Be — ?■ 'K^iKi decays are sensitive to the value of 9k, which 
will be tested by the running LHC and forthcoming SuperB experiments; 

5. The LPFs is larger than 80% for almost all decay modes. That means that these 
pure annihilation decays of Be meson are dominated by the longitudinal polarization 
fraction. 

These charmless hadronic Be meson decays will provide an important platform for 
studying the mechanism of annihilation contributions, understanding the helicity struc- 
ture of these considered channels and the content of the axial-vector mesons. 
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Appendix A: Wave functions and distribution amplitudes 



For the wave function of the heavy Be meson, we adopt the form (see Ref. [10|, and 
references therein) as follows, 

^bA^) = ^ [(f + mBj750B.(x)],^ . (Al) 

where the distribution amplitude (pBc is of the form [36] in the nonrelativistic limit, 

M^) = ^5ix-me/mBj . (A2) 

In fact, we know little about (pBc for heavy Be meson. Because of embracing b and c 
quarks simultaneously. Be meson can be approximated as a non-relativistic bound state. 
At the non-relativistic limit, the leading 2-particle distribution amplitude (pBc can be 
approximated by delta function [36], fixing the light-cone momenta of the quarks according 
to their masses. According to Ref. [36], this form will become a smooth function after 
considering the evolution effect from relativistic gluon exchange. 
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For the wave function of axial-vector meson, the longitudinal (L) and transverse(T) 
polarizations are involved, and can be written as. 



^li^) = -^i^{mAiA4>\{x) + 



(A3) 



{x) + mAie^upal^l^e*r^nPv''(t)\{x)]^^ , (A4) 



L T 

where e^' denotes the longitudinal and transverse polarization vectors of axial-vector 
meson, satisfying P ■ e = in each polarization, x denotes the momentum fraction carried 
by quark in the meson, and n = (1, 0, Oy) and v = (0, 1, O-r) are dimensionless light-like 
unit vectors. We here adopt the convention e°^^^ = 1 for the Levi-Civita tensor e^'^°'^. 

The twist-2 distribution amplitudes (Pa{x) and Eqs. (lA3|A4p can be parame- 

terized as (23. 129I: 



= - x) 

0^(2;) = ^x{l-x) 



%A + 3aiA (2a; - 1) + 4^ l^i^i'^^ - 1)^ - 1) 
4^ + 3ai^ {2x - 1) + 4^ ^(5(2x - 1)^ - 1) 



(A5) 
(A6) 



Here, the definition of these distribution amplitudes (f)A{x) and satisfy the following 

normalization relations: 



1 

0ip^(x) = a[Jip^ 



2^6' 







2^6 ■ 



(AT) 



where ajjsp^ = 1 and agiPi ~ -'- have been used. 

As for the twist-3 distribution amplitudes in Eqs.f lA3|A4p . we use the following 
form [29] : 



<P'Aix) = ^ l^aU^x -ir + l a^A (2x - l)(3(2x - 1)^ - 1) 
0aW = {^(1 - 3^)(«i4 + a^Ai.'^x - 1))} . 



(A8) 
(A9) 



0a(^) = ^ + (2x - 1)^) + aL(2x - 1)^ ^ 



3/a 

4^6 6^3; 



x(l - x)(aQ^ + a1^(2x - 1)) 



(AlO) 
(All) 



where Ja is the decay constant of the relevant axial-vector meson. When the axial-vector 
mesons are Kia and Kib-, x in the distribution amplitudes stands for the momentum 
fraction carrying by the s quark. 
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The Gegenbauer moments have been studied extensively in the hteratures (see Ref. 24 
and references therein), here we adopt the following values: 



"2/1 



''2/8 



"IK, 



'■2K, 



-0.02 ±0.02 
-0.04 ±0.03 
-0.07 ±0.04 
0.00 ±0.26; 



'1/1 



a 



1/8 



-1.08 ±0.48; a! 



OKiB 



0.02 ±0.10; 



a 



IKiB 



-1.04 ±0.34; 
-1.06 ±0.36; 
-1.11 ±0.31; 
-0.05 ±0.03; 
= 0.14 ±0.15; 
0.17 ±0.22. 



"Ihi 



'1/18 



-1.95 ±0.35; 

-2.00 ±0.35; 

-1.95 ±0.35; 
= 0.08 ±0.09; 
= -1.95 ±0.45; 



(A12) 
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